Transport behaviour of a Bose Einstein condensate in a bichromatic optical lattice 
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The Bloch and dipole oscillations of a Bose Einstein condensate (BEC) in an optical superlattice 
is investigated. We show that the effective mass increases in an optical superlattice, which leads to 
localization of the BEC, in accordance with recent experimental observations [16]. In addition, we 
j^Q 1 find that the secondary optical lattice is a useful additional tool to manipulate the dynamics of the 

' atoms. 
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I. INTRODUCTION 



Interference pattern of intersecting laser beams create a periodic potential for atoms, which is known as an optical 
lattice [l| . Ultracold bosons trapped in such periodic potentials have been widely used recently as a model system for 
the study of some fundamental concepts of quantum physics. Josephson effects _|2|], squeezed states [3j, Landau-Zener 
tunneling and Bloch oscillations [1] and superfluid-Mott insulator transition [flare some examples. An important 
promising application under study is quantum computation in optical lattices [6]. Optical lattices are therefore, of 
' particular interest from the perspective of both fundamental quantum physics and its connection to applications. 
| Using superposition of optical lattices with different periods [3], it is now possible to generate periodic potentials 
characterized by a richer spatial modulation, the so-called optical superlattices. The light-shifted potential of the 
y | superlattice is described as 

^) = ^lCOS 2 ^)+y 2 COS 2 ^ + ^. (1) 

i Here d\ and di are, respectively, the primary and secondary lattice constants, V% and V2 are the respective amplitudes 

and 4> is the phase of the secondary lattice. When <fi — 0, all sites of the lattice are perfectly equivalent due to 
the symmetries of the system, so that the population and onsite energies are same at each site. An asymmetry is 
k*" ' introduced when <fi ^ and hence the onsite energies are not the same at each site. 

Theoretical interest in optical superlattice started only recently. These include work on fractional filling Mott 
insulator domains ||, dark Q and gap [l(| solitons, the Mott-Peierls transition [ll[, non-mean field effects [Hj] 
and phase diagram of BEC in two color superlattices [HI . In a recent work, the analogue of the optical branch in 
solid-state physics was predicted in an optical superlattice [l4j]. Rousseau et al. [l5| have considered the effect of a 
t^. ' secondary lattice on an one dimensional hard core bosons (strongly correlated regime). A detailed theoretical study 
of the Bloch and Bogoliubov spectrum of a BEC in a one-dimensional optical superlattice has been done In a 
very recent experiment [13], it was observed that the center of mass motion of a BEC is blocked in a quasi-periodic 
lattice. Considering the fact that these optical superlattices are now being realized experimentally and interesting 
experiments are being done routinely, we were motivated to study the influence of the secondary lattice on Bloch 
oscillations and dipole oscillations of atoms. 



II. BLOCH OSCILLATIONS 



We consider an elongated cigar shaped BEC confined in a harmonic trap potential of the form V/j (r, z) = ^ (io 2 r 2 - 



lu 2 z 2 ) and a one-dimensional tilted optical superlattice of the form V op (z) — Er (si cos 2 (^) + S2 cos 2 (||)) +mgz. We 

have taken a particular case of di — 2d\ — 2d. Here s± and S2 are the dimensionless amplitudes of the primary and the 

- - ' h 2 ir 2 Er 

secondary superlattice potentials with s\ > S2- Er — is the recoil energy ( lor = — — is the corresponding recoil 

2md 2 a 
frequency) of the primary lattice. We take lo t » uj z so that an elongate cigar shaped BEC is formed. The harmonic 

oscillator frequency corresponding to small motion about the minima of the optical superlattice is u s w — — — . The 

peak densities in each well match the Gaussian profile. Since the array is tilted, the atoms undergo coherent Bloch 
oscillations driven by the interwell gravitational potential mgz. The BEC is initially loaded into the primary lattice 
and the secondary lattice is switched on slowly. The frequency of each minima of the primary lattice is not perturbed 
significantly by the addition of the secondary lattice. Here u s >> lu z so that the optical lattice dominates over the 
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harmonic potential along the z-direction and hence the harmonic potential is neglected. The strong laser intensity will 
give rise to an array of several quasi-two-dimensional pancake shaped condensates. Because of the quantum tunneling, 
the overlap between the wave functions of two consecutive layers can be sufficient to ensure full coherence. We study 
now the Bloch dynamics of the BEC in the tilted optical superlattice by solving the discrete nonlinear schroedinger 
equation (DNLSE). The dynamics of the BEC is governed by the Gross-Pitaevskii equation (GPE), 



lh it = ~ 2^ v ' c + { y,io(r ' z) + Vop{z) + 3olc|2 } c 



(2) 



where <?o 



Anh 2 a 



with a the two body scattering length and m the atomic mass. In the tight binding approximation 



the condensate order parameter can be written as 

CM) = 



3 



(3) 



where Nt is the total number of atoms and <f>(r — rj) — <f>j is the condensate wavefunction localized in the trap j with 
J dr(j)j<j)j +1 ss 0, and J dr \(j)j\ 2 = 1; ^fj(t) is the j th amplitude. *ffj(t) = pj (t) exp (i6j (t) ) where pj = with Nj 
and 6j being the number of particles and phases in the trap j respectively. Substituting the Ansatz ([3]) in ([2|) , we find 
that the GPE reduces to the DNLSE, 



89 a 
'~dt 



= ~\ {(1 - a{-iy- 1 ) *i- x + (1 - a(-iy) + ( £j + A |^-| 2 ) 



(4) 



Here e j = j-J dr [«L (V^) 2 + (V ho (r) + V op (z)) \^\ 
~J dr \&V<l>j ■ Wj+i + <f>j {Vho(r) + V op (z)) <j> i+i 



Ji 2 = Jn ± % where Jn w 



\-2 



exp 



A = So^r J dr l^l 4 , a = One can show using Jj = 

that there are distinctly two Josephson coupling parameters, 
and A ■^■s 2 exp ^— ^ [161] . We have rescaled 



time as t 



dU, Ej = wbJ, where lob 



mgXi 



is the frequency of Bloch oscillation and Ai is the 



wavelength of the laser creating the primary lattice. In order to understand the Bloch and dipole oscillations, we solve 
the DNLSE using a variational approach adopted from [l8[. The Hamiltonian function corresponding to the DNLSE 
Eq. (g|) reads 



H 



E 



i 



{(i-(-i)V) 



3+1 



) + (i - (-lr^) + + ej |*/ + | i^i 4 



(5) 

where J2j = 1- To analyze the Bloch dynamics, we study the dynamical evolution of a site dependent Gaussian 
wavepacket, which we parameterize as 



*,•(*) = VKexp 



(6) 



where £(t) and j(t) are, respectively, the center and width of the condensate, p(t) and 5(t) are their associated 
momenta, and K (7, £) a normalization factor. Here (— is the phase of the wave packet at the j th site. Clearly, 
depending upon whether j is odd or even, the phase is ±?. As explained in ref.(15), as the condensate moves from 
one well to the next, it acquires additional phase, which depends on the height of the barrier. As the height of the 
barrier alternates, the phase also alternates. 

The dynamics of the wave packet can be obtained by the variational principle from the Lagrangian, L = y\ i^fj^J* — 

H, with the equations of motion for the variational parameters = £, j,p, 8, <p given by ^f^ 1 = |pS The phase is 
used to enforce a constraint. The Lagrangian is derived as 

A 



L 



• 7 2 <5 
"8" 



20F7 



+ {cos 4>cosp + a sin (j) sin p} exp(— 77) — V(j, £), 



(7) 



where V = ^ + and V( 7 , £) = K dj e 3 exp (-2^) 
The variational equations of motion are derived as: 



dV 



(8a) 
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Figure 1: Oscillations of the center of mass £(t) is depicted for two different values of the secondary lattice strength, a — 0.1 
and a = 0.4. The other parameters are £(0) = 0, p(0) = 0, 5(0) = 0, 7(0) = 10, A = 20, lub = 2. On increasing the strength of 
the secondary lattice, the amplitude of the center of mass motion reduces. 



£ = [cos^sinp — a sin cos p] exp(— 77), 



5= [cos (f>cosp + a sin ifisnip] exp(— rj) 



2A 



4dV 



7 3 Y / 7r 7 97 ' 
7 = 7<5 [cos(/>cosp + a sin <fi sin p] exp(— 77), 



(8b) 

(8c) 
(8d) 



tan0 = a tan p. 

Since cos 2 <j> + sin 2 tj> = 1, together with equation (|8aH8ep . we get the following constraints on cos</> and sin^: 



cos ( 



cosp 



V 



cos 2 p 



o ^2 ' 

or sin p 



a sinp 



cos 2 p - 



a 2 sin 2 p 



(8e) 

(9a) 
(9b) 



Corresponding to the variational equations (|8aj|8ep and constraints (|9afl9b[) the effective Hamiltonian is written as 



H 



A 



20F7 



— cospy 1 + a 2 tan 2 pexp(— 77) + V(j, £). 



(10) 



We first study the Bloch oscillations. For the tilted periodic potential the on-site energies are written as £j = jujb- 

Using equations (|8a,fl8e|) . we find V = £lub and p — —lob- We solve the variational equations of motion numerically 
for the following initial values £(0) = 0, p(0) — 0, 6(0) — 0, 7(0) = 10 and the parameters A = 20, ojb = 2. The result 
for the center of mass £(i) is depicted in figure 1 for two different values of the secondary lattice strength, a = 0.1 
and a = 0.4. Clearly on increasing the strength of the secondary lattice from a — 0.1 to a = 0.4, the amplitude of 
the center of mass motion reduces. The secondary lattice serves to break the discrete translational invariance of the 
system, thus favouring localization of the wave function. Optical superlattices with higher periodicities will block the 
center of mass more strongly. The observed damping (with respect to time) in fig. 1 is due to interactions. In the 
absence of interactions, the center of the BEC for po — goes roughly as £(£) —(1 — a 2 )(l — cosuist), while in the 

presence of interactions, the oscillations roughly decreases as w —(1 — a 2 ) (l — cxp ^— ^pr^j cosajB^j. Here, 7/ 

is some final value of 7. Clearly when there is no interaction, there is no damping of the Bloch oscillations in time 
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Figure 2: A plot between m* and a shows that as the strength of the secondary lattice increases, the effective mass also 
increases. Therefore the origin of the localization of a BEC in an optical superlattice is due to an increase in the effective mass. 

but there is a reduction in the amplitude by a factor (1 — a 2 ) due to the presence of the secondary lattice. In order 
to understand the origin of this blocking of the center of mass motion, we derive the effective mass (to*)" 1 = as, 

(l + a 2 tan 2 p) 3/2 cxp(r7) 

m (i 2 1 4 Wi 5T- ( U ) 

cosp (1 — or tan^pj (1 — cr) 

A diverging effective mass m* — > oo as t — > oo due to interactions leads to a self-trapping of the wave packet [17]. 
In the expression for the effective mass (eqn.ll), in the absence of interaction, the factor exp(7y) is constant since 7 
tends to a final value 7/ and S(t) w So (initial value). This can be seen from equations 8c and 8d. The effective mass 
is now enhanced due to the presence of the secondary lattice. Since A = 0, the effective mass stays constant in time 
and the Bloch oscillations show reduced oscillations compared to the case for a single frequency optical lattice but 
does not show damping in time. On the other hand when A ^ 0, and t — ► 00, 7 — > 7^ and S(t) w -|y= i so that 

to* — > 00. This causes not only a reduction in amplitude but also damping in time. It is interesting to note that, we 
now have an additional handle to tune the effective mass. A plot between to* and a (for p = 0) in fig. 2 shows that as 
the strength of the secondary lattice increases, the effective mass also increases. Therefore the origin of the reduction 
of the amplitude of Bloch oscillations of a BEC in an optical superlattice is due to an increase of the effective mass. 
Dynamics of localized excitations, such as solitons depends on the effective mass, hence the secondary lattice emerges 
as a useful additional handle to manipulate localized excitations. 




III. DIPOLE OSCILLATIONS 

We study now the dipolc oscillations. Instead of the gravitational potential, we consider a sufficiently large (u) z w u> s ) 
magnetic harmonic potential superimposed on the optical lattice, Sj = £lj 2 , where ft = —f Q — • The variational 
equations of motion give V(7,£) = £1 {^j- + and p = — 2Q£. In the regime of negligible mean field interaction 

(A = 0) and small momenta p, the equation for the center of mass is £(t) = (l — en 2 ) p. Consequently, the center 
of mass obeys the equation of an undamped harmonic oscillator, £ = where the frequency of dipole oscillation, 
uj 2 , = 20 (l — a 2 ) = lo 2 (^t) is reduced in the presence of the secondary lattice since to* > to. We consider 
the initial conditions £(0) = and p(0) = p - The center of mass in the A = regime and small momenta is 

(l-a 2 ) 1/2 

£(f) w i — -j=L — sinujdt- In the low momenta limit, the amplitude of the center of mass decreases with increasing 

1/4 

. In the experiment of ref. [16], ut z w 2ir x 10 
Hz and Ai w 830 x 10~ 9 nm. This corresponds to a very low value of £1 w 0.0001 (in dimensionless units). 
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Figure 3: A plot of the dipole oscillations for a = 0.1 and a = 0.7. The other parameters are £(0) = 0, p(0) = 0.1, 5(0) = 
0, 7(0) = 40, A = 5, Q, = .0002. We notice that increasing the strength of the secondary lattice, the dipole oscillations are 
blocked, in accordance with the experimental observations of [l7|. Since, we are in the negligible mean field interaction regime, 
the dipole oscillations are not damped. 



We solve the variational equations of motion numerically for the following initial values: £(0) = 0, p(0) = 0.1, 6(0) = 
0, 7(0) = 40 and the parameters: A = 5, fi = .0002. The result for the dipole oscillation is depicted in figure 3 for 
two different values of the secondary lattice strength a = 0.1 and a = 0.7. For A = 5, we are still in the regime of 
negligible mean field interaction and we do not expect any damping. On increasing the strength of the secondary 
lattice, the amplitude of the center of mass £(t) is reduced in accordance with the experiments of fl7| . This reduction 
in the amplitude of the dipole oscillation on increasing the strength of the secondary lattice is due to an increase in the 
effective mass, as mentioned earlier in this paper. The initial value of the effective mass can be positive (cospo > 0) or 
negative (cospo < 0). Let us suppose that cospo > and initial values: 7(0) = 70, 6(0) = 6q = and £(0) = Co = 0. 



The initial value of the Hamiltonian is Hn = 



77-^= cosn 

2V7T70 ^ 



\Jl + a 2 tan 2 po exp (- 



^j 2 - . Since the Hamiltonian 



is conserved, we have Ho — Tr-r* — 
helps to keep H > 0, therefore, 



- cospo\/l + a 2 tan 2 po exp {^\l 2 — - L ~i~^J + ^4"- The parabolic external potential 



A fi 7 2 

27^ + ^-^ >0 - 



The trajectories in the 7 — 6 plane are given by 



(12) 



87 2 log 



.Oil 



-H 



cos po yj 1+ck 2 tan 2 pg 



r 



(13) 



Fig. 4 shows a plot of the center of mass for a = 0.1 , f (0) = 0, p = 0.1, <5 = 0.1, 70 = 10, A = 47, 57, uj = 0.0002. 
We notice that for such high values of A, the dipole oscillations are completely blocked. Both interactions and 
secondary lattice induced disorder cooperate to block the center of mass motion. For the higher value of A, the center 
of mass stops at an earlier time, which again is in accordance with experiments [l~7| . From equation (13), we notice 
that 5 — > 00 as t — » 00. Therefore, for large time, 



t « (1 - a 2 ) sinp exp ( - %^ ) -> ( 14 ! 
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Figure 4: Center of mass motion for a = 0.1 , £(0) = 0, p = 0.1, <5 = 0.1, 70 = 10, A = 47, 57, w = 0.0002. We notice that for 
such high values of A, the dipole oscillations are completely blocked. Both interactions and secondary lattice induced disorder 
cooperate to block the center of mass motion. 



The center of the BEC wavepacket stops and the effective mass goes to infinity and there is an energy transfer 
from the kinetic energy to the internal modes, since 6 is the momentum associated with the width 7. This is the self 
trapped regime. We also find that the final value of center of mass £/ is not the same as £o- For a fixed A, an increase 
in the secondary lattice potential will block the center of mass at an earlier time. 



In conclusion, we have studied the Bloch and dipole oscillations of a Bose Einstein condensate trapped in an optical 
superlattice. In particular, we have shown that due to the addition of the secondary lattice, the center of mass motion 
is blocked which leads to a blockage of the center of mass motion. This effect is due to an increase in the effective 
mass in the presence of the secondary lattice. The frequency of the dipole oscillations is also found to be reduced due 
to the secondary lattice. These results are in accordance with recent experiments 17]. The secondary lattice is found 
to be a promising tool to investigate and manipulate localized excitations. 
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IV. CONCLUSIONS 
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